Large N expansion for the 2D Dyson gas 



A. Zabrodin * P. Wiegmann ^ 

December 2005 



Abstract 

We discuss the 1/N expansion of the free energy of A'^ logarithmicahy interacting 
charges in the plane in an external field. For some particular values of the inverse 
temperature /3 this system is equivalent to the eigenvalue version of certain random 
matrix models, where it is referred to as the "Dyson gas" of eigenvalues. To find the 
free energy at large A'^ and the structure of l/A'^-corrections, we first use the effective 
action approach and then confirm the results by solving the loop equation. The 
results obtained give some new representations of the mathematical objects related 
to the Dirichlet boundary value problem, complex analysis and spectral geometry 
of exterior domains. They also suggest interesting links with bosonic field theory 
on Riemann surfaces, gravitational anomalies and topological field theories. 
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1 Introduction 



In this paper we discuss the large expansion of the A^-fold integral 

Zm= \Ar,{z,)ril^^'^^''^d'z,. (1.1) 
i=i 

Here is a function of z = x + iy and z = x — iy, (3 > is a. parameter, A]\f{zi) = 
Yii^jizi — Zj) is the Vandermonde determinant, and (Pz = dxdy. The "Planck constant" 
U is introduced here to stress the quasiclassical nature of the large N limit: N ^ oo 
together with — so that to = is kept fixed. 

The integral is equal to the partition function of the statistical ensemble of 2D 
Coulomb charges in the external potential W (the "Dyson gas" PP). Different aspects of 
2D Coulomb plasma were discussed in j^-IZ], for further applications see 0-^21 ■ In this 
interpretation, the parameter (3 is the inverse temperature. At large 13 the Dyson gas is 
believed to form a Wigner crystal. In this paper we assume, however, that (5 is such that 
the system is in the liquid phase as A^ — > oo. At some particular values of (3 the model 
can be viewed as the eigenvalue version of certain ensembles of random matrices [TH] 
(normal or complex matrices at /? = 1 jElE] and normal self-dual matrices at /? = 2). 

An important information is encoded in the 1/A^-expansion of the model. For random 
matrices, it has the meaning of the "genus expansion" since the 1/A^ order of perturbation 
theory graphs is determined by their Euler characteristics. In the sequel, we prefer to 
work with the equivalent /i-expansion, thus emphasizing its semiclassical nature: 

logZ^ = c{N) + § + + Fi + OiU) (1.2) 

The explicit form of the c{N) is given below. It can be absorbed by normalization. With 
some abuse of terminology we call fi^ log a free energy. 

When A^ becomes large new macroscopic structures emerge. The gas segregates into 
"phases" with zero and non-zero density separated by a very narrow interface. Let 
D be the domain in the complex plane where the density is non-zero (it may consist 
of several disconnected components). In the first approximation, the gas looks like a 
continuous charged fluid trapped in the domain D. The density at any point outside it 
is exponentially small as A^ — oo. 

The first two terms in (jl.2|) . Fq and -F1/2, are of purely classical nature in the sense 
that only the static equilibrium state of the charges (the saddle point of the integral) 
contributes to them. The leading contribution to the free energy, Fq, is basically the 
Coulomb energy of the charged fluid in the domain D. Taking into account the discrete 
"atomic" structure of the Dyson gas, which implies a short- distance cutoff and entropy 
of macroscopic states, one is able to find the correction F1/2 to the "classical" free en- 
ergy. The next term, Fi, apart from further corrections of the classical nature, includes 
contribution from small fluctuations about the equilibrium state. 

The (3 = 1 Dyson gas confined to the line is related to the model of Hermitian random 
matrices. The 1/A^-expansion of this model beyond the leading order has been obtained 
in the seminal paper ^Hl- Recently, there was a progress in understanding these results 
from the algebro-geometric point of view ^[ 02] and in extending them to other 
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matrix models |2m 1^ 1^ . In the genus-1 correction was interpreted in terms of 
free bosons on Riemann surfaces. Our results for Fi (partially reported in [221) ^njoy a 
similar interpretation and lead to interesting connections with spectral geometry of planar 
domains. In particular, our results suggest a formula for the determinant of the Laplace 
operator in exterior planar domains in terms of the conformal map of the domain onto 
the exterior of the unit circle. For polynomial potentials, Fi enjoys a finite determinant 
representation similar to the one known in topological field theories. 

Our results for Fi suggest a new deep connection between the 2D Dyson gas and 2D 
quantum gravity. This connection does not explore the well known approach to random 
surfaces through a scaling limit of random matrices. It rather indicates that the density 
of 2D Dyson particles can be treated as a fluctuating 2D metric. We do not develop this 
approach in this paper. 

In the rest of the introductory section we fix the notation and present some standard 
exact relations to be used in the sequel. We follow [22 | l2 H l25]. 

The main observables in the Dyson gas statistical ensemble are mean values and 
correlators of symmetric functions of the particles coordinates. Let A{zi, . . . , zj^i) be 
such a function, then the mean value (A) is defined by 



1 r ^ 



i=i 

A particularly important example is the density 

p{z) = hY^5{z~z,) (1.3) 

3 

where 6{z) is the two dimensional 5-function. Instead of correlations of density it is often 
convenient to consider correlations of the field 

V{z) = -P j\og\z-C?p{Od\ (1.4) 

from which the correlations of density can be found by means of the relation 

ATi(3p{z) = -Aip{z) (1.5) 

Here and below, A = Adzdg is the Laplace operator. Clearly, (p is the 2D Coulomb 
potential created by the charges. 

Handling with multipoint correlation functions, it is customary to pass to their con- 
nected parts. For example, in the case of 2-point functions, the connected correlation 
function is defined as 

{p{Zi)p{z2))^ = {p{Zi)p{z2)) - {p{Zi)) {p{z2)) 

The following variational formulas hold true: 

^P^'^^ = ' ' ^P^'^^P^'^^^'^ = ' mW = 6Wiz,)6Wiz2) ^'-'^ 

These formulas are exact for any finite A^. They follow from the fact that variation of 
the partition function over a general potential W inserts X)j ^{^^ ~ ^i) into the integral. 
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More generally, the connected part of the {n + l)-point density correlation function is 
given by the linear response of the n-point one to a small variation of the potential. 

Let f = f{z, z) be a function in the complex plane (for brevity we write simply f{z) in 
what follows). Summing over the charges, we get the symmetric function f{zi) = tr/, 
where the notation is inspired by related models of random matrices. Mean values and 
correlators of such functions are expressed through those of densities: 

(tr/) = j {p{z))f{z)d'z, (tr/itr/2), = J {p{z{)p{z,)) J,{z,)Uz2)d' z,<fz, (1.7) 

and so on. 



2 Ward identities 

Clearly, the integral remains the same if we change the integration variables. In 
other words, it is invariant under reparametrizations of the 2;-coordinate. This leads to a 
number of Ward identities which are our basic tool to compute the free energy at large 
N. 



2.1 Holomorphic form of the Ward identity at finite N: loop 
equation 

We begin with a holomorphic reparametrization Zi ^ Zi + e{zi). Let us apply it to the 
integral (jl.lj) . = J e~R^(^i'---'^'^) flj (Pzj, where the energy is 

- n-'E = pJ2hg\z, -z,\ + n-'Y.^i^j) (2-1) 

In the first order the integrand transforms as 

fdE , , dE—^ 

while the volume element Jlj (^^j undergoes the scaling (Weyl) transformation 



-^3 
3 



\^Y.^de{z,)^dt{z{] 
I 



■^3 
3 



The invariance of the integral is then expressed by the identity 

d 



valid for any e. Introducing a suitable cutoff at infinity, if necessary, one sees that the 
2D integral over Zi can be transformed, by virtue of the Green theorem, into a contour 
integral around infinity and so it does vanish. 

Let us take 

t{z,) = ^— , (2.2) 

Z - Zi 
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where z is a complex parameter. The singularity at the point z does not destroy the above 
identity since its contribution is proportional to the vanishing integral § dzij {zi — z) over 
a small contour encircling z. We explore this singularity in the sequel. Therefore, we 
have the equality 



E 



-^^ + 7 



Z-Zi \Z-ZiY 



where Tl ^d^^E = —f3 ^ h ^dW{zi). Using the identity 



l-^i Zi Zi 



1 2 1 

^ {z - Zi){z - Zj) ^. {z - Zi){zi - Zj) i{z- 



we rewrite it in the form 

(T) = 0, (2.3) 
where we define the holomorphic component of the stress energy tensor 



z - Zii^ 



(2.4) 



This is the holomorphic form of the Ward identity. In terms of the field '^[^z^ p.4|l it 
The correlator at coinciding points is understood as ({d(p{z)Y) = lim {d(p{z) d(p{z')). 



We have got an exact relation between one- and two-point correlation functions, valid 
for any finite A^. For historical reasons, it is called the loop equation. Since correlation 
functions are variational derivatives of the free energy, the loop equation is an implicit 
functional relation for the free energy. However, it is not a closed relation. It can be 
made closed by some additional assumptions or approximations. A combination with the 
1 /A^-expansion is particularly meaningful. 



2.2 Path integral representation of the partition function and 
Weyl form of the Ward identity 

Here we give somewhat heuristic but transparent arguments to calculate the free energy 
by using the path integral representation and the Weyl form of the Ward identity. The 
results will be justified and refined by means of the loop equation. 

Density as a metric. At large A^, when the density can be treated as a smooth 
function, the partition function of the Dyson gas can be represented as a path integral 
over densities p{z). Symbolically, we write the partition function as 

Zn = J [Dp]e- M^f^ (2.6) 

where the action A is to be determined. In this approach p{z) dzdz appears as a metric 
gabdx"'dx^ written in the conformal gauge Qah = Sabp{z), while the loop equation has the 
meaning of the Ward identity with respect to holomorphic diffeomorphisms. 
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Stress energy tensor. The action in Eq. ()2.6|) can be determined by its response to a 
change of the metric. The response of the action to variation of the metric is generated 
by the stress energy tensor (s.e.t.). Under the Weyl transformation p —>■ p + 6p, the 
change of the action is 

-6A=- I bp p-^ (2.7) 

TT J 

where T^z is the trace of the s.e.t. The easiest way to determine T^j is to use the 
conservation law of the s.e.t. reflecting the reparametrization invariance of the action. 
In the conformal gauge the conservation law reads 

37 ^pd{p-^%,)=^ (2.8) 

where T is the holomorphic component of the s.e.t. It has been already derived directly 
from the finite dimensional integral (ll.lj) . The result (Eq. ()2.4|) ) can be rewritten as 

-2/5r(z) = 2/?|^^p(C)rf\+(a(^(z))2 + (2-/3);.c?V(^) (2.9) 

Let us apply d to this equality. Taking into account that diXj z) = 7t6{z) and using p.5|l . 
we obtain dT = tt p (j)(p — dW + |(2 — (3)d\og pj . The conservation law then states that 



(2.10) 



where A is the integration constant which will be determined by the normalization con- 
dition / (p) cPz = hN. 



The action. Integrating Eq. (j2.7|) with T^z given by (j2.10|) . one determines the action^ 
up to a constant: 

-A = ^f (f3-^<f{p-^A)if + SnW - 47r{2-f3)h\og{e^p)) p£z (2.11) 

Here p^^A is the invariant Laplace-Beltrami operator in the metric pdzdz. 

Let us discuss the physical meaning of this action. It consists of energy and entropy 
contributions. The energy is given by (j23])- On the scales much larger than the mean 
distance between the charges the system can be treated as a charged liquid with the 
electrostatic energy 

- hEo[p\ = P j j p{z) log \z - CI p{z')d^zd^z' + j W{z)p{z)dh (2.12) 

= ^J {r'v{p^'A)^ + 87rW)pd'z. 

It gives the leading term of the action. The correction |(2 — /5) / plogpd'^z results from 
the discrete "atomic" structure of the Dyson gas. The argument below goes back to 
Dyson pOj. 

^The effective action for some other models of random matrices has been discussed in Refrs. |26ll27j . 
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The subleading term of the action consists of two contributions of different nature. 
One correction comes from the sum I]j^jlog|zj — Zj\, when passing to the continuous 
theory. Namely, one should exclude the terms with i = j, writing 

J2^0g\z^-Zj\ =J2\og\Zi- Zj\ -J2\ogi{zj) 

or 

'^\og\ Zi — Zj\ = fi^"^ j p{z) \og\z — z'\p{z')d'^zcfz' — J p{z) log i{z) d"^ z 

where £ is a short-distance cutoff (which may depend on the point Zj). It is natural to 
take the cutoff to be 

. /j: (2.13) 

which is the mean distance between the charges around the point z. This gives the 
improved estimate for the electrostatic energy: 

E[p] = Eo[p] -^pj plogpd'z + ^phNlogh- j,hN (2.14) 

where 71 is a numerical constant which can not be determined by this argument. 

Another correction comes from the integration measure when one passes from the 
integration over Zj to the integration over macroscopic densities. We can write 

Y[d'zj = N\ J[p] [Dp] 
j 

where [Dp] is an integration measure in the space of densities, J[p] is the Jacobian 
of this change of variables and the factor A^! takes into account the symmetry under 
permutations (all the states that differ by a permutation of the charges are identical). 
To estimate the Jacobian, we divide the plane into microscopic "cells" such that j-th 
particle occupies a cell of size i{zj), where i{zj) is the mean distance ()2.13|) between the 
particles around the point zj. All the microscopic states in which the particles remain in 
their cells are macroscopically indistinguishable. Given a macroscopic density p, J[p] is 
then approximately equal to the integral JceiisHj ^"^^jy with each particle being confined 
to its own cell. Therefore, J[p] ~ Ilji'^{zj), and thus log J[p] (sometimes referred to as 
entropy of the state with the macroscopic density p) is given by 

h\ogJ[p] = - J pXogpd^z + hN\ogh + -f2hN (2.15) 

where 72 is another numerical constant^. 

The subleading term of the action, — | (2-/3) / p log pd'^z, is thus the sum of the 
contribution due to the short distance cutoff and the entropy contribution. They cancel 
each other in the ensemble of normal self-dual matrices (at (3 = 2). 

^Combining (|2.14(l . I|2.15|l and taking into account the factor iV! in the measure, we obtain for the 
c{N) in Eq. H1.2|l : c{N) = log A^! + ^(2 — /?) log h + jN, where 7 = 71 + 72 is a numerical constant. 
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The Weyl form of the loop equation. The gravitational Ward identities state that 
the expectation value of the variation of the action vanishes. They can be written in two 
complimentary (holomorphic and Weyl) forms 

(T) = 0, (T,,) = 0. (2.16) 

The first identity, generated by the holomorphic component of the s.e.t., is the holomor- 
phic loop equation ()2.5p . The second one, — ^ {%z) = (pf^) — ; 

{^p)-W (p) + \{2-P)h (plog(e^+V)) = (2.17) 

is another form of the loop equation. 

The Weyl form of the Ward identity represents the gravitational anomaly. It is based 
on the fact that the diffeomorphism ()2.2|) is not holomorphic at the points Zi (positions 
of the particles). At these points de{z) = e7T6{z — Zi) and the holomorphic component T 
of the s.e.t. is not analytic. It has simple and double poles. 



3 The structure of the expansion 

The large N limit we are interested in {Nh = to remains finite) corresponds to a very 
low effective temperature of the gas, when fluctuations around equilibrium positions of 
the charges are negligible. The main contribution to the partition function then comes 
from a conflguration, where the charges are "frozen" at their equilibrium positions and, 
moreover, the gas can be treated as a continuous fluid at static equilibrium. Mathemati- 
cally, all this means that the integral (jl.ip is evaluated by the saddle point method, with 
only the leading contribution being taken into account. Fluctuations around the saddle 
point give 1/A^ corrections. 

The path integral ()2.6|) makes the structure of the 1/A^ expansion intuitively clear. 
First we flnd an equilibrium ("classical") density pd which minimizes the action A[p]. 
Then, separating the classical part of the density, p = Pd + ^Sp we can write 

A[p] = A[pd] + y / 5piz) Kiz, z')Spiz')d^z + ... 

where 

K{z, z' 



bp{z)bp{z') 



P=Pcl 



is the kernel of an integral operator K. The path integral representation of the Dyson 
gas immediately produces the flrst two leading contributions to the free energy: 

h^logZ^ = -A[pd] - -logdeiK +0{n'') (3.1) 

The flrst term is the classical value of the action, the second one is due to the Gaussian 
fluctuations around the classical solution'^. 



■^One can obtain this results by a direct iteration of the Ward identity (|2.17|l without appeaHng to 
the path integral representation. 
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Since the leading part of the free energy is the classical value of the action, the 
connected part of the pair density correlation computed in the leading order is equal to 
the kernel of the operator inverse to K: 

h-' {p{z)p{z'))^ = {5p{z) 5p{z')) = k~\z, z') in 0) (3.2) 

The term coming from fluctuations (the second term in 1)3.11) ) is not the only contribution 
to Fi. The action itself and the solution Lp^i of the equation 

W-^,i = -{2-(3)h\og ( -^A^d j (3.3) 

minimizing the action depend on Ti. Expanding the density, pd = Po + + • • •, we 
obtain 

A[pcl] = A[po] j Pl/2{Z)K{Z, Z')py2iz')£zd^z' + ..., (3.4) 

The action evaluated at po yields the first two leading terms of the free energy. Up to a 
constant we write 

-A[po\ =Fo + hFy2 ■ 
The second term in fl3.4|) contributes to the next order: 

-A[pd]+A[po] = n'F}'\ 



Summing up, we obtain three leading orders of the expansion ()1.2)1 : 



P Poiz)log 



PoiOd'zd'C 



Fl/2 



-i(2 



j3) / pologporf z 



F, = ) + F| 



(0) 



--log detK, 



(1) 



Pi/2{z)K{z, z')pi/2{z')d'^z(fz' 



(3.5) 

(3.6) 
(3.7) 
(3.8) 

(3.9) 



They are expressed entirely through the mean density and the pair correlation function 
()3.2|) computed in the leading order. This remains to be the case for the higher order 
corrections as well. 

The mean density and the pair correlation function were found in [21]. Using these 
results one is able to compute the free energy by means of ()3.5|) - ()3.9p . We summarize 
the results in the next section. In the sequel we derive them on a more solid basis by a 
direct iteration of the holomorphic loop equation. 

We note also that the ^-expansion of the free energy can be written in the "topo- 
logical" form Y,g>o^'^^ Fg, where each term has its own expansion in e = (2 — P)h: 
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4 The main results 



We introduce the following notation: 

aiz) = -^AW{z), x{z) = log^/^), a = ^-^ (4.1) 

4.1 Summary of the results 

In order to compute the first three leading contributions to the free energy we need the 
following results. 

• The mean density computed as a power expansion in h vanishes outside a 
bounded domain D (we assume that D is connected). Inside the domain the first 
two orders are 

po = /3-V, zeD; (4.2) 

py2 = -{2-(3)K-\ zeD. (4.3) 

The leading correction to the density pi/2 is singular at the boundary (see Eq. (|5.16|) ). 
Corrections exponential in h make the density a smooth function falling exponen- 
tially outside D. 

• The shape of D is determined by the potential W and the number of particles, 
as is described below (see fj4.16p ). It is the subject of the inverse potential problem 
in 2D. 

• The pair correlation function of densities is related to the Dirichlet boundary 
problem in the complimentary domain C \ D (the exterior of D). Given a function 
f{z), let f^{z) be its harmonic continuation from the boundary to the exterior of D 
(the solution of the exterior Dirichlet boundary value problem). The pair correla- 
tion function of densities (and thus the kernel K{zi, Z2)) is completely characterized 
by the integral relation 

47r/5 / f{z)k~\z,z')giz')d'zdh' = - j^^^V{f - nV{g-g'')d'z + j^^f^gd'z 

^\ (4.4) 
for any smooth functions /, g. Hereafter, K ^{z,z') is the kernel of the operator 
K~^. The first integral goes over the exterior of the domain D while the second 
one is over the entire plane. Note that / — vanishes on the boundary. 

• The spectral determinant det K in ()3.8|) is a ratio of two spectral determinants of 
the invariant Laplace-Beltrami operators e~'^^A in the conformal metric vro" = e^^: 

det(-e-2xA(p,p) 
logdeti^^log . (4.5) 

One of them acts on the entire plane (and so does not depend on D). The other one 
is the Laplace-Beltrami operator in the exterior domain C \ D acting on functions 
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vanishing at the boundary. The determinants of the Laplacians in exterior, un- 
bounded domains were introduced in The definition is more involved than the 
usual zeta-function regularization. For the details see |2Hj and references therein. 
However, explicit formulas for the exterior determinants were not known. Our re- 
sult for Fi'^^ obtained by a direct solution of the iterated loop equation suggests 
such a formula (see below). 

To present the results for the free energy we need some more notation. Let ds = \dz\ 
be the line element along the boundary of D, k be the local curvature of the boundary, 
and e'^dwdw = dzdz be the conformal metric induced by the conformal map w{z) of the 
exterior of D onto the exterior of the unit circle. 



Free energy: 



2-P 



'1/2 



2/3 



1_1 a{Odhd\ 
z C 

cr log(7rcr) d'^z 



(4.6) 
(4.7) 



(0) 



247r 



/ \V{(f) + x)\'d'w-2<f {(f) + x)\dw\ 

J\w\>l J\w\=l 
ilT JC OTT JdD 



^ 47r 



C\D 



247r 



(4. 



|Vx|V^ 



+ 



27r JdD 



dnXds (4.9) 



Note that F1/2 here differs from ()3.6|) by a constant times that is a matter of 
normalization. We have determined the coefficient /x only at /? = 1. In this case 
12 = -1/8. 

Special cases: 



- A quasiharmonic potential: W = —\z\'^ + V{z) + V{z). The density is uniform 
inside the domain, a = l/vr, F1/2 ~ ^^^^ 

1 



Fr = 0, 

{(l)dn(l) + 2(f))\dw\ 



(4.10) 



247r J\w\=i 

-(3 = 2. In this case the "classical" corrections (-F1/2 and -Fi^'*) vanish for any 
potential W . 

pio) 

is given by ()4.8|) . 

These results, being combined with the spectral determinant representation ()3.8|) . 
fl4.5|) . suggest an explicit formula for the determinant of the Laplace-Beltrami operator 
in exterior domains. For the interior domains, such a formula has been known due to 
|29| l30j. It is called the Polyakov-Alvarez formula. For simply-connected domains, in can 
be written in terms of the conformal map of the domain onto the unit disk. (A sketch of 
the derivation is given in Appendix D.) Our results suggest that for the exterior case this 
formula remains basically the same, with the interior conformal map being substituted 
by the exterior one. Below we list some formulas of this type. Their rigorous proof is a 
challenging problem which is beyond the scope of this paper. 
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The Polyakov- Alvarez formula for exterior determinants (a conjecture): 

det(— e'^^Afp, „) 1 /■ 1 r 1 r 

det(— A(j-<^p) 127r jC\d ott J9d 47r jC\c 



(4.11) 

(the divergent terms proportional to the area and perimeter are omitted). In the 
l.h.s., det(— A(Q^j-,) is the regularized spectral determinant of the Laplace operator 
in the Euclidean metric. This determinant is expressed entirely in terms of the 
boundary value of the metric e'^dwdw = dzdz induced by the conformal map w{z) 
of the exterior of D onto the exterior of the unit circle: 

logdet(-Ac\D) = T^/ 0(/t + e-*)rfs = -^/ (09„0 + 20) . (4.12) 

12tT JdD ll-K J\w\=l 

Eqs. (imi) . KT^ give: 



logdet(-e-2^Ac\D) = " / \V{<P + x)? d^w - 2 I {<P + x)\dw 

IZTT [J|«)|>1 J\w\=l 

Determinant formula (a conjecture). As a by-product we suggest another 
interesting formula for the spectral determinant. Let tk = ^qq z'^'zdz be har- 
monic moments of the exterior of the domain D and nto be its area. Let us assume 
that the domain is such that all the moments at k > m vanish. Then, up to a 
constant, the exterior spectral determinant is 

where Fq is given by ()4.6|) with a = l/vr. Although the r.h.s. looks like a complex 
number, it is actually real, as it will be clear from the derivation below (see ()(i.7|l . 



4.2 Comments and details 

The mean density in the leading order. In the leading order Eq. reads 
2/3 J log \z - C\poiC)d\ + Wiz) = A , A = ^(2 - f3)hX 

It must be valid in the domain where p 7^ 0. We assume that this domain (the support 
of the density) is bounded, otherwise the normalization condition, i.e., that the integral 
/ pd'^z be finite, can hardly be satisfied. We denote the support of the density by D. 
In this paper, in order to avoid additional technical complications, we consider only the 
case of connected domains D. 

Upon taking the 2;-derivative of the equation above, we get, for z inside the domain 

D: 

difo{z) = dW{z) {z e D) (4.14) 
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with 

^o(^) = -/5/log|^-ClVo(C)A 

being the 2D Coulomb potential created by the equilibrium ( "classical" ) configuration of 
the charges characterized by the density po- This equation just states that the total force 
experienced by a charge at any point 2;, where po{z) 7^ 0, is zero. Indeed, interaction 
with the other charges, d(pQ{z), is compensated by the force dW{z) due to the external 
field. 

The solution is conveniently expressed through the function a ()4.H1 . For the model 
to be well-defined, we assume that a{z) > and tends to a positive constant as \z\ — >• 00. 
Applying dz to both sides of eq. ()4.14j) . one obtains the solution: po{z) = a{z)/[3 inside 
D and Pq{z) = outside it. In terms of the potential ipo this solution reads 

M^) = - [ \og\z-C\MOd\ (4.15) 

Assuming, without loss of generality, that G D and W{0) = 0, we fix A = —(po{0), and 
so W{z) = ipo{z) —ipo{0). Plugging this into ()H.fi|l . we find the leading contributions 

to the free energy ()4.6|) . ()4.7p . Fq is the electrostatic energy of the domain D charged 
with the density a{z) and with a point-like compensating charge at the origin. Different 
aspects and applications of the functional ()4.6p were discussed in jH], [211 

The most nontrivial part of the problem is to find the shape of D. It is determined 
by eq. ()4.14|1 and by the normalization condition. Using the Cauchy integral formula, 
we can write these conditions in the form 

JdD Z — C 

(4.16) 

^ JD 

To solve them for D provided W{z) and to are given, amounts to a version of the inverse 
potential problem in two dimensions. We assume that the potential W is such that the 
solution exists and is unique. In general, the solution is not available in an explicit form. 
In this paper we do not address this question. Our goal is to express corrections to the 
free energy ()4.6|) in terms of the domain D. 



Pair correlation functions. The correlation functions in the leading order can be 
found using the general variational formulas ()1.6|) . where the exact free energy is replaced 
by Fo: 

Basically, these are linear response relations used in the Coulomb gas theory [^j. In 
this approximation, the 2D Coulomb plasma is represented as a continuous charged 
fiuid, so the information about its discrete microscopic structure is lost. This is a good 
approximation at distances much larger than the mean distance between the charges. 

Here are the results for the correlation functions of the type (jl.7p . The mean value 
is obvious from the result for the mean density given above: 

n(3(tTf)= I a{z)fiz)(fz +0{h) (4.17) 

JD 
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The variation w.r.t. the potential yields the connected parts of pair correlators 
To present the result, we need some elements of the Dirichlet boundary value problem. 
Given a function f{z), let f^{z) be its harmonic continuation from the boundary of D 
to its exterior, i.e., the solution of the exterior Dirichlet boundary value problem. The 
solution is given by the formula 

f''{z) = -^{ f{0dnG{z,0m (4.18) 
Zn JdD 

where G{z,^) is the Green function of the domain C \ D. It is the symmetric function of 
two points uniquely determined by the following properties: 

A,G{z,C) = 27r6{z-C) in C \ D , ^(^,0=0 if z e dD 

For simply connected exterior domains C\D the Green function can be expressed through 
the conformal map w{z) from C \ D onto the exterior of the unit circle^: 



G'(^,C)=log 



w{z) - w{C) 



1 — w{z)w{() 



(4.19) 



As C ^ z, it has the logarithmic singularity G{zX) log k ~ C|- The connected pair 
correlator ()1.7|) is then given by 



A'K(3{iifiig)^ = A'K(3 [ f{z)k-\z,z')g{z')d^zd?z' = f VfVgd^z-i fd^g'^ds (4.20) 

J Jd JdD 

where V/ is the gradient of the function /. Alternatively, with the help of the Green 
formula Eq. ()4.2()|1 can be rewritten as a bulk integral ()4.22j) . where the integration goes 
over the exterior domain and the entire plane: 

A7T/3 J f{z)k-\z,z')g{z')d^zd''z' = - j^^^V{f-f'')V{g-g")d^z + J^VfVgdh 

(4.21) 

In particular, for the connected correlation functions of the fields ^p{zi), (p{z2) this formula 
gives (if zi^2 e C \ D): 

{^iz^Mz2)), = Giz,, Z2) - G(zi, oo) - G(oo, z^) - log + O^U) (4.22) 



2/3^ 

where logr = lim2_»oo(log \z\ + Gi^z^ oo)) is the (external) conformal radius of D 



Spectral determinants. The relation ()4.5|1 between the spectral determinants, 

log det k = log det(-e-2^Ac^p) - log det(-e-2^Ac) , (4.23) 

can be understood using formula (j4.2(jp . The presence of the bulk and boundary terms 
in the r.h.s. suggests to separate boundary and bulk values of the functions / and g. 
Specifically, we write f = fn + f , where fn is the harmonic continuation of the / to the 

^Throughout the paper, the map w{z) is normahzed as w{z) = z/r + 0{1) as z ^ cx) with a real r. 
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interior of D while the boundary value of / is zero, and similarly for g. Using the Green 
theorem, one can see that the bulk and boundary contributions completely separate^: 

4nf3 J f k-^g = jj (-e^^^^A) g e^'^dh + £^ fNg ds 

Here 

Ngiz) = d:^gH{z) - d-g'^iz) 

is the Neumann jump operator [d^ are normal derivatives from inside and outside respec- 
tively) which sends a function on the boundary to the difference of normal derivatives 
of its harmonic continuations inside and outside the domain D. This means that the 
operator is the direct sum of the Laplace operator in D (with the Dirichlet b.c.) and 
the Neumann jump operator. Therefore, det{K^^) ~ det(— e^^-^Ao) det iV, and so 

= i log det(-e2>^AD) + ^ log det N (4.24) 

(we omit an irrelevant factor). This expression can be simplified by means of the following 
relation between the properly regularized functional determinants pHj : 

logdet(-e"^^AD) + logdetiV + logdet(-e"2^Ac^p) = logP(D) + logdet(-e"^^A(c) 

Here -P(D) = ^/¥o ds is the perimeter of D in the metric vrcr. 

This relation, rigorously proven in |2H] , is known in the mathematical literature as the 
"surgery formula". It is clearly motivated by the "cut and paste" physical arguments. 
Consider the free bosonic theory in the whole plane with the quadratic action 5*0 = 
/ iVXpd^z (for brevity, we do not indicate the metric explicitly). The path integral 
l\DX\ exp(— S'o[X]) is equal to (det(— A(j-^))^^/^. On the other hand, let us fix a domain 
D and represent the action as 

Sq = I \VX?d^z^ I IVXIV^ 
Jd 7C\d ' ' 

Decompose the field X inside D into the sum of the field X such that X = on c)D and 
the harmonic field Xu'- X = X + Xh- Let X = X + X^ be the similar decomposition 
for the field X outside D (we then have Xh = X^ on dD). It is easy to see that these 
fields separate in the action as follows: 

So[X] = - I XAXd^z- ( XAXd^z+ <f XuNXuds 
Jd JC\d JdD 

This separation implies the surgery formula. The term logP(D) is due to the zero mode 
of the Neumann jump operator which we did not take into account. 

Hence, Fi'\ the "quantum" part of Fi, is to be identified with | log det(— e~^'^A(Q) — 
I logdet(— e~^^A(j-^^p), where the determinants have to be properly regularized. Some 
additional efforts are required to refine these arguments. In particular, one should justify 
the choice of the background metric in D and take care of the zero mode of the Neumann 
jump operator. This is beyond the scope of the present paper. The next section provides 
an alternative derivation of Fi. 

^We deliberately keep the factor e^^ — ira in some formulas in order to emphasize the fact that the 
Laplace-Beltrami operator is an invariant operator with respect to the metric ttct written in the Weyl 
gauge. Although the metric cancels in the formula below, it appears in the spectral determinants. 
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5 Corrections to the free energy from the loop equa- 
tion 



In the previous sections we have found that the asymptotic expansion of the partition 
function as — > has the form 

= N\h^^'-^^^ e-^^ exp + ^ + + 5: h^-'F,/^ (5.1) 

The corresponding expansions of the mean values of p and are 

{p{z)) = po{z) + hp^/2{z) + T^p^{z) + O^K") (5.2) 

{^{z)) = <^o(^) + h^i/2{z) + n\^{z) + 0{h^) (5.3) 

where pi{z) = SFi/6W{z) and (pi{z) = -(3 j Pi{C) ^og\z - Cpc?^C- The terms Fq and F1/2 
are given by ()4.6|1 and ()4.7|1 respectively. The higher corrections are due to fluctuations 
of the charged particles around the equilibrium state. In principle, they can be found by 
expanding the loop equation ()2.5j) 

^ / ?}^LmM01 = ((a^(,)).) + (9V(=)) (5.4) 

in powers of Tl and solving the inhomogeneous linear integral equations obtained in this 
way. A similar approach has been developed in the case of Hermitian matrix ensembles 
[in]. This is what we are going to do in this section. We restrict ourselves by F1/2 and 
Fi. Calculations of higher order corrections are rather tedious. 



5.1 Iteration of the loop equation 

As it was already pointed out, the main contribution to the partition function as ^ — > 
comes from a configuration, where the charges are "frozen" at their equilibrium positions. 
Correspondingly, the averages take their "classical" values {^p{z)) = fo{z), and multipoint 
correlators factorize in the leading order: {dip{z) d(p{z')) = dipo{z)dipo{z'). Under this 
assumption, the loop equation becomes a closed relation for ipo: 

^J?^^^^^C-{0M4 (5.5) 

where we have omitted the last term in (j5.4j) which is of the next order in h. Let us 
apply dz to both sides of the equation. This yields: dW{z)AipQ{z) = dipo{z)A{pQ{z). 
Since AipQ^z) cx po{z) we obtain 

po{z) [dMz) - dW{z)] = (5.6) 

This equation should be solved with the additional constraints / po{z)(Pz = to (normal- 
ization) and po{z) > (positivity). The equation tells us that either dipo^z) = dW{z) or 
Po{z) = 0. Applying d^, we get Aip^^z) = AW{z). This gives the solution for po and ipo 
already obtained in Sec. 2 by less formal arguments (see (j4.15|) ). 
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Now we are in a position to develop the /i-expansion of the loop equation ()2.5|) . First 
of all, we rewrite it identically in the form 

^ I Liz, (A^(C)) d\ = idM^)r - {d {{^{z)) - Mz) 



which is ready for the ^-expansion. Here 

,(..C)^ ^'^K)-^^°M (5.7, 

is the kernel of the integral operator in the l.h.s. (the "loop operator"). The zeroth order 
in h gives equation (|5.6|) which implies the familiar result ifo{z) = — /□ log \z — (\'^a{()d'^( 
for the ipo- To proceed, one should insert the series (|5.3p into the loop equation and 
separate terms of order h, etc. The terms of order h and h'^ give: 



1 
1 

2^ 



L{zX)^^i/2{Qd\ = -(2-/5)aVo(^) 

(5. 



L{zX)^Vi{Qd\ 



dvi/2{z)] +(2-/3)aVi/2(^) 



where 



uj{z) = lim 



u;{z) 

(5.9) 



h ^ lim {d(p{z) d(p{z')) 

is the connected part of the pair correlator at merging points. 

The expansion of the loop equation can be continued order by order. In principle, 
this gives a recurrence procedure to determine the coefficients (pk{z). However, each step 
requires solving integral equations in the plane, that is not easy to do explicitly. Another 
difficulty is that in general one can not extend these equations to the interior of the 
support of the density because the ^-expansion may break down or change its form there. 
Indeed, in the domain where the density is macroscopically nonzero, the microscopic 
structure of the gas becomes essential, and one needs to know correlation functions at 
small scales. Nevertheless, at least in the first two orders in h the equations above can 
be solved assuming that z G C \ D. Note that in this region all the functions (pk{z) are 
harmonic. If these functions are known, the corresponding expansion coefficients of the 
free energy in can be obtained by "integration" of the variational formulas ()1.6|) . 



5.2 The solution for F1/2 

We start with the order h. We need to solve the first equation in ()5.8|) . Using ()1.5|) . 
()1.6|) . we rewrite it in the form 

Lia,z)d^z='-^ ipl^ (5.10) 



SW{z) ' ' ' 2p JD{a-z 

where a is an arbitrary point in C \ D. Using the variational technique developed in 
(see also Appendix A), one can verify that F1/2 given in eq. ()4.7|) does solve this 
equation. 
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Fi/2 from the loop equation. Here we give some details of this calculation. Exactly 
the same scheme is used in the next subsection while solving the loop equation for Fi. 



It is convenient to use the notation ()4.1|1 x(^) = log JTia{z) in terms of which F1/2 



— Jo e^^xd'^z- The variation of this functional reads 



SF 



1/2 



Here 



2-P 

Tip 



Sn{z 



JdD JD 

dn{SW{z)-6W"{z)) 
4:TTa{z) 



(5.11) 



(5.12) 



is the normal displacement of the boundary (Fig. 1) under variation of the potential, 
with the convention that 6n > for outward displacement (for the proof see [211 1^ ) 
and 

A6W{z) 

^X{z) = -- — ^ 5.13 
87rcr(z) 

Note that the l.h.s. of ()5.10p has the meaning of the variation of F1/2 under a small change 
of the potential proportional to L{a, z), 6W{z) oc L{a, z) (a plays the role of a parameter). 
Therefore, the result is given by (j5.11|) where 5n and 5x are taken from (j5.12j) . (j5.13p with 
5W{z) = L{a,z). It remains to plug the explicit form of the L{a,z) ()5.7|) and simplify 
the result. In the course of this calculation, the frequently used formulas are: 



A,L(a, z) = Aird. 



a{z) 



(5.14) 



and 



9D 



fiz)dJL{a, z)-L''(a, z)) ds = 2tii I ZifMfl 

JdD a — z 



(5.15) 



(for any smooth function /). It is implied that z G D, a G C \ D. In these formulas, the 
Laplace operator and the harmonic continuation are applied to z. We have: 



6W{z] 



L(a, z) d z 



2-P 
47r/3 



x(z)dJL(a,z)-L^(a,z))ds - - / (2x(z) + l)AL(a, z)d^z 

dD 2 JD 



2-£ 



X{z)cr{z) 



dD Z 



dz+ / i2xiz) + m 



cr(z 



d'z 



After transforming the first (contour) integral to the integral over the domain D, 



X[Z)(T[Z} 



dz = —2i / d. 



dD Z 



a 



'x{z)a{zy 



dh 



one can see that the result is indeed equal to the r.h.s. of (|5.1(J|) . 
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Figure 1: The normal displacement of the boundary. 



The results for pi/2 and ^1/2- The corrections to the mean values of p and (p are 
conveniently expressed through the function x ()4.H) and its harmonic continuation x^i^) 
from the boundary of D to its exterior. In terms of these functions 



6W{z) Att/S 



0(z; D)Ax{z)-6{z;dD)dM^) 



'X''iz))-l5'{z;dD) 



(5.16) 

Here Q{z; D) is the characteristic function of the domain D (1 if 2; G D and other- 
wise), 6{z; dD) is the 5-function with the support on the boundary (/ f{z)6{z] dD)d'^z = 
§Q^f{z)\dz\ for any smooth function /) and 6'{z;dD) is its "normal derivative", i.e. a 
function such that / f{z)6'{z; dD)(Pz = — /^q dnf{z)\dz\. The singular function pi/2 is to 
be understood as being integrated with any smooth test function. The first term in the 
r.h.s. describes the change of density in the bulk. The second one (the "simple layer") 
describes a shift of the boundary of the domain D. The last one, the "double layer", 
means smoothing off the edge of the density support (for j3 ^ 2). Using the operator K 
introduced in Sec. 3, the correction pi/2 can be written as 

p„2{z) = -{2-m-\]{z) - ^<5'(z;9D) 

At /3 = 2 the correction pi/2 vanishes. The first correction to {f{z)) reads: 



-(2-/5) [x{z) - x"" (oc) + 
-i2-(3) [x^(z) - x^(oo)^ 



z eD 

zeC\D 



(5.17) 



Due to the double layer, this function is discontinuous across the boundary. Therefore, 
for z e C \ D we have: ^(^1/2(2) = -{2- l3)dx^ (z). 



5.3 The solution for Fi 

Now we have to solve the second equation in ()5.8p . where dipi/2 = —{2—[3)dx^ and uj{z) 
is still to be found. If the point z is in C \ D, then eq. ()4.22|) yields j^ : 

{d^{z) d^{z'))^ = 2(3h''dA' {G{z, z') - log |z - z'\) + OiK")) 

Since the r.h.s. is regular for all z, z' G C \ D, the points can be merged without any 
difficulty and the result does not depend on the particular limit z' — ^ z. Using the 
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expression of the Green function ()4.19|) through the conformal map w{z), we obtain 



12 lim d.dJGiz, z') - log \z -z'\)= ^ - | f ^ ) = {w, z] 

where we use the standard notation for the Schwarzian derivative. The function lo is thus 
given by 

uj[z) = ^{w; z} (5.18) 
After plugging the above result for ipi/2, the second equation of ()5.8p acquires the form: 

/ Wt) "^'^ = l(9x''{z)? - dh^'iz)] + ^{w; z] (5.19) 

where a = — \J~^ (as is in (j4.ip ). We consider this equation for 2; G C \ D, where the 
both sides are harmonic functions. 

Our strategy is as follows. We make a guess for Fi and then show that the variational 
derivative obeys ()5.19p . More precisely, let P"^"^ be a trial functional of D (and thus 
of W) represented as an integral over the domain D or its boundary. The form of the 
trial functional is suggested by the path integral arguments of the previous section. 
Since Fi has dimension 0, we consider only dimensionless functionals. We want to find 
/ dF^^^"^ / 5W {z)L{a, z) (P z and compare with the r.h.s. of ()5.19p (again, a is a point in 
C \ D). The latter quantity can be computed by the method outhned in the previous 
subsection. 

A comment on the meaning of the calculations below is in order. They consist of 
recognizing that the l.h.s. of eq. (j5.19p is a variation of Fi over the holomorphic compo- 
nent of the metric and then applying the Polyakov- Alvarez formula (|4.1ip . We plan to 
elaborate on this point elsewhere. 

The structure of eq. (j5.19j) suggests to find different terms of the solution separately. 
Let us start with the term proportional to a^. Here are the main steps of the calculations. 
Consider the functional 

I^'^ = ^ (/^ Y^X?d'z - £ xdnX'^ds^ (5.20) 
Its variation is (see Appendix B): 

5/« = -^/ \d^{^-x")\'5nds + ^i 6xdn{x-x'')ds~^ [ 6xAxd'z (5.21) 
An JdD L -I 27r Jao 2n Jd 

In the same way as for F1/2 one can check that 



2 



6W{z] 



L{a,z)d''z={dx''{a)f-d\''{a) 



so Fx is expected to contain the term a^/^^-* (compare with given by (j4.9|l ). This 
contribution is of the "classical" nature since no fluctuations of the charges positions are 
taken into account. 
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Next, we consider 

^^'^ = -^f i'Pdn<P + 2(t))\dw\ (5.22) 

ITI J|«,| = l 

(we remind that 0(ui) = log |-2'(w)| where z{w) is the conformal map from the exterior of 
the unit circle onto C \ D inverse to the w{z)). The variation of this functional is found 
in Appendix B (eq. 



5/(2) = ^j> (v\z){w- z} + u^iz){w; z} - 2k\z)) 5n{z) ds (5.23) 

Here 

v{z) = \w'iz)\^ (5.24) 
is the normal unit vector to the boundary® and 

w{z) 



k{z) = dn log 



is the local curvature of the boundary (counted w.r.t. the outward pointing normal 
vector). Using the rules explained above, we get: 



L(a, z) a z = : f az 



6W{z) Am JdD a — z 

Now consider the functional 



1^'^ = ^ I \^X?d'z + -<l Kxds (5.25) 



with the variation 



51^^^ = ^<f> (\Vx?-2dlx+2^^d^x)5nds+- I (n + d^xWds - - j AxSxdh 

ZtT JdD ^ ^ Vr JdD TT Jd 

The variational derivative SI^^^SW looks rather complicated. However, for its convolu- 
tion with L{a, z) one obtains a surprisingly simple result 

5/(^) . / N ,2 1 /■ f^ds 1 f + iK' 



L{a, z) d z = f = — : f dz 



6W{z) ' 27r JdD {a — zY 27ri JdD a — z 

(note that x cancels!). Combining it with the corresponding result for J(2), we obtain: 



5(/(2)_j(3)) I . v^{w-z} + v^{w-z} + 2ik' 

L{a,z)d^z = : i> — ! ! dz 



5W{z) ' Alii JdD a — z 

1 / i^'^{z){w;z} 1 r {w;z} 

dz = — : f dz = {w[a); a} 



2^1 JdD a — z 2T[i JdD a — z 

where fl8.2p has been used. 



^It is worthwhile to mention here the useful formulas for normal and tangential derivatives: 9„/ 

dz 



^dzf + vdzf , dsf = ivdzf — ivdzf and for the line element on the boundary curve: ds = -£fr^ = iv{z)dz. 
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At last, the functional 



with the variation 



2tt Jd 



(5.26) 



277 JdD 

is a "zero mode" of the loop operator: 

5/(0) 



6W{z] 



L{a, z) d z = 



Summing all the contributions with appropriate coefficients that follow from the r.h.s. 
of the loop equation, we thus find Fi. 

F, = ) + = + 1(J(2)_ j(3)) + ^j(o) 

where the coefficient /i can not be determined from the loop equation restricted to C \ 
D. (In the whole plane, the loop operator does not have zero modes but we have to 
be restricted to C \ D because of unknown properties of correlation functions at small 
distances in the bulk.) For (3=1 the coefficient fi can be fixed by comparison with the 
explicit solution for centrosymmetric potentials (see Appendix C): 

= --^ (at /3 = 1) 
o 

With the help of the Green theorem, we present the result in the form appearing in 
SecEl 

^ ' \V{(P + x)\'d^w-2f {(l) + x)\dw\ 

\w\=l 



Fi 



2A7T 



w\>l 



+ — 

4:71 



I |Vx|V^- / xdnX"ds 
Jd JdD 



(5.27) 



where x the first two terms is regarded as a function of w via x = x{^{f^))- Let us 
also list some equivalent forms: 

1 



Fi 



2Att J\w\=i 
1 -6a2 
247r 



(09„0 + 20)|rf^/;| 

/ \Vx?d^z + 2i Kxds 
Jd JdD 



(5.28) 



a 



Fi 



1 -6^2 
247r 



\w\>l 



f X{dnx" + 2K)ds + ^ i dnXds 

JdD 2tT JdD 

\V{(P + x)\^d^w-2i i(P + x)\dw\ 



w\=l 



C 



a 

+ — 

Ait 



V(0 + x)|Vz 

/ \V{<P + x'')\'d'w-2i {<p + x'')\dw\ +^fAxd'z. 

J\w\>l J\w\=l /TT JD 



(5.29) 
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6 Models with quasiharmonic potentials 



The models with W of the form 

W = -\z\'^ + V{z) + V{z) 

("quasiharmonic potentials") generalize the Ginibre-Girko ensemble ^H]- Note that for 
potentials of this form the integral diverges unless V is quadratic or logarithmic 
with suitable coefficients. The simplest way to give sense to the integral when it diverges 
at infinity is to introduce a cut-off, i.e., integrate over a suitably chosen big but finite 
domain in the plane. Then the large expansion is well-defined and the results for the 
general potential presented above still make sense. For details and rigorous proofs see 

m- 

In the case of quasiharmonic potentials the formula for Fi drastically simplifies since 
the function x vanishes, and so only the first integral in ()5.28|) survives: 

^1 = - ^ / (^^-^ + 20) \M (6.1) 
ZAtt J\w\=i 

In the particular case /3 = 1, W{z) = —zz the formula yields = log to that 
coincides with the result of j3] obtained by a direct calculation. 

According to the conjecture of Section 4.1, eq. ()6.H) can be understood as the formula 
for the regularized determinant of the Laplace operator A(j~<^j^ = Ad^dz in the exterior 
domain C \ D with the Dirichlet boundary conditions: 

Fi = -ilogdet(-AcJ (6.2) 

The first term is the bulk contribution (for the metric induced by the conformal map 
it reduces to a boundary integral), while the second term is a net boundary term. The 
"classical" contribution F^^^ to Fi vanishes in this case. 



6.1 The case of rational dV{z) (quadrature domains) 

There is a special class of domains for which our result (j6.H) can be made more explicit. 
Consider domains such that z'(w) is a rational function, 



z'{w) = r Y[ 



w - a,- 



i=o 

In the mathematical literature, they are called quadrature domains |36i- One can show 
that quadrature domains are density supports for the models with potentials such that 
dV{z) is a rational function. All the points a, and bi must be inside the unit circle, 
otherwise the map z{w) is not conformal. As w oo, z{w) can be represented as 
a Laurent series of the form z{w) = rw + uq + 0{w~^). On the unit circle we have 
\dw\ = ^ and 4>{w) = ^(\ogz'{w) + log z' {w~^)) , where the first and the second term 
(the Schwarz reflection) are analytic outside and inside it, respectively. (Recall that 
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z{w) = z{w) and our notation z'{w ^) means dz{u)/du at the point u = w Plugging 
this into we get: 



247ri 



f log z'(w 

J\w\=l 

If , __^dw 



]-du,\og z\w) + — 
2 w 



dw — 



, logz'iw ,^^} dw 
w 4:8m J\w\=i z'[w) 



The integrals can be calculated by taking residues either outside or inside the unit circle. 
The poles are at oo, at 0, and at the points and 6j. The result is 



^ 2'(ai)=0 z'{bi)=oo 



(6.3) 



If the potential V{z) is polynomial, V{z) = J2T=i'^kZ^, i-e., = as > m for some 
m > 0, then the series for the conformal map z{w) truncates: z{w) = rw + Z^S^ UiW~^ 
and 



m— 1 



z{w) = r W_{1 



aiW 



j=0 



is a polynomial m w ^ (all poles 6j of z'{w) merge at the origin). Then the last sum in 
()6.3|) becomes m log r and the formula ()6.3|) gives 



i^i = -^iog|r^ n 



m— 1 



2'{aj)=0 



(6.4) 



This formula is essentially identical to the genus-1 correction to the free energy of the 
Hermitian 2-matrix model with a polynomial potential computed in 



6.2 The determinant formula. 



For polynomial V{z) Fi enjoys an interesting determinant representation. Set 



Dr 



det 



dtodtjdtk 



0<j,k<m-l 



where Fq is the leading contribution to the free energy regarded as a function of to and the 
coefficients tk (and the complex conjugate coefficients tk). We need the residue formula 
for the third order derivatives of Fq |37j : 



d^F, 



hj{w)hk{w)hi{w) dw 
w 



dtjdtkdti 2711 J\w\=i z'{w)z'{w~^) 
Here hj{w) is the following polynomial in w of degree j: 

hj{w) = (-2"'(^))+ for J > 1 and hQ{w) = 1 



(6.5) 
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where (...)+ is the positive degree part of the Laurent series. Using this formula, we 
compute: 

^ _ 1 / dwo f dwm-i det [hj{wj)hk{wj)] 

{2TCi)'^ J\wo\=l Wo "' J\w^.i\=l Wm-l UZ~0^ Z'{wi)z'{w^^) 

Clearly, the determinant in the numerator can be substituted by ^det'^{hj{wk)) and 
det [hj{wk)] = (m — 1)! r5™('"~^)Am(wj), where Am{wi) is the Vandermonde determinant. 
Each integral in ()6.6|) is given by the sum of residues at the points Oj inside the unit 
circle (the residues at = vanish). Computing the residues and summing over all 
permutations of the points Oj, we get: 

It is not difficult to see that Hi^i = (— l)™m(m — l)r™~^tm (we regard tm, the last 
nonzero coefficient of the V{z), as a fixed parameter) and so we represent Fi ()6.4|) in the 
form ^ ^ ^ 

Fi = — log Dm - — (m^-3m + 3) logr - — (m-1) logt„ + const (6.8) 

where const is a numerical constant. Since O^^Fq = 2 logr (see |211I2S]), we obtain that 
Fi, for models with polynomial potentials of degree m, is expressed through derivatives 
of Fo. 

According to the conjecture of Section 4.1, the formula below suggests a new rep- 
resentation of the spectral determinant of the Laplace operator. Up to a constant we 
have 

log det (—A(r^n) = — log det ( -r-^r—^ — I (m^— 3m+3) „° (m— 1) logt^ 

2 ^ V 24 ^ ydtodtjdtk J 24^ ' dtl 24^ ; & m 

(6.9) 

where j, run from to m — 1. Similar determinant formulas are known for genus- 1 
corrections to free energy in topological field theories ^B]- However, they have not been 
identified with spectral determinants. 



7 Spectral determinant of the Laplace-Beltrami op- 
erator for the Dirichlet problem 

To give an interpretation of the result for Fi, we recall the formula for the spectral 
determinant of the Laplace-Beltrami operator in a compact bounded domain M in the 
plane (assumed to have topology of a disk). For the derivation see j2n[Eni, Section 1 of 
PU] and Appendix D. The Laplace operator acts on functions vanishing on the boundary 
(the Dirichlet boundary conditions). Being written in the Weyl gauge, Qah = e^*5afe, it 
has the form 

where w is a holomorphic coordinate on M. The Polyakov-Alvarez formula gives the dif- 
ference between the spectral determinants of the Laplace operators in the metric e^'^dwdw 
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and in some fixed reference metric. Assuming tliat tlie latter is just tlie standard fiat 
metric dwdw in tfie plane, the formula reads (see (4.42) in the first paper in [^): 



logdet(-e-2*AM) = -7^/ e^'^d'w + -^<f e*|rfz/;| + ^ loge 

4:7ie^ Jm ^\/T^t JdM 6 



Aire 
1 

127r Jm 



(7.i; 



67r JdM 



k^ldwl — I A^d'^'w 



An Jm 



Here e is an ultraviolet cutoff and k is the curvature of the boundary w.r.t. the reference 
metric. The first three terms diverge as e — 0. A sketch of the derivation is given in 
Appendix D. 

We are going to show that our result for the "quantum" part of the free energy, Fi^\ 
given by ()4.8|) . ()5.27|) - ()5.29|) agrees with ()7.1|) generalized to exterior domains. Indeed, 
let us adopt this formula for the exterior of the domain D. For this purpose, we map it to 
the exterior of the unit circle Uext by the conformal map z ^ w{z) and choose M in 1)7.11) 
to be the exterior of the unit circle, M = Uext, with the metric $(w) = xi^i^)) + (piw). 
Then 

e^^dwdw = e^'^dzdz (7.2) 
Taking into account that for the exterior of the unit circle k = — 1, we rewrite ()7.H) as 

logdet(-e-'<Ac,„) = -ij /^^^ ad^. + i £ + 1 log. 

(^\dw\ 



127r 



V(^\'^d^w - 2 



\w\>l 



w\=l 



Att J\w\>1 



(7.3) 

On the complex plane C (without boundary) with coordinate z in the metric e"^^ we have: 
logdet(-e-2^Ac) = -^ + Uoge - ^ J^lVxl'd'z (7.4) 



Therefore, we can write: 

det(-e"^^A(j-) 



log 



^0 ^ 1 , P(D) 
— + -loge--^ 
4e^ 6 Ay/ne 



+ 



12-K 



det(-e-2xA(^^J 
/ |V$prf2«;-2/ ^\dw\-?,i dnx\dw\ 

J\w\>l J\w\=l J\w\=l 



(7.5) 



1 



12tx JC 



\Vx\'d'z 



The first two divergent terms can be absorbed by normalization and so we ignore them 
in what follows. The third one is proportional to the perimeter of the density support 
(in the metric e'^^): 

P(D) = i y^ds 

JdD 

It can be directly verified that the perimeter functional obeys the relation 

6P{D) 



6W{z] 



L(a, z) d z = 



{aeC\D) 
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so all the three divergent terms are "zero modes" of the loop operator (for the first two 
this is obvious). 

Comparing ()7.5|1 with ()5.27|1 . we can represent F\ in the form 



1 det'(-e-2xA(j-) 



^'-2'°^det'(-e-xAc;„;-^ 



47r 



dD 



+ 



fx 



2tt JdD 



dnXds 



(7.6) 

Here logdet' means the finite (as e ^ 0) part of ()7.5|) . It is the "quantum" part of the 
answer. Obviously, Fi is related to gravitational anomalies. This relation awaits further 
clarification. 



8 Appendices 

A. Some useful formulas 

Here we fix the notation and collect some formulas often used in the main text. 



Integral formulas. Throughout the paper, all contours are assumed to be anticlock- 
wise oriented and the normal vector looks outward D. 

• Cauchy's integral formula (/ is any smooth function): 

f fiQdC _ 1 f d^fiQd'C ^ i -fiz) , ZED 
2m JdD z — ( TT Jd z — 

• The Green formula: / fAgd'^z = — / Vf^gd^z^ f fdng\dz\. 

Jd Jd Jod 

• The Hadamard variational formula: the variation of the Green function of C \ D 
under small deformation of the domain D with the normal displacement 6n{z) is 

5Giz,, Z,) = ^i drM^u OdnGiz2, OHO M^l 

Zn JdD 

Variation of contour integrals. Consider the contour integral of the general form 

/ F{f{z),dnf{z))ds 

JdD 

where F is any fixed function. Calculating the linear response to the deformation of the 
contour (described by the normal displacement 6n{z)), one should vary all items in the 
integral independently and add the results. There are four elements to be varied: the 
support of the integral / , the 9„, the line element ds and the function /. By variation 
of the f we mean integration of the old function over the new contour. This gives 
§ 5ndnF ds. The change of the slope of the normal vector results in 5{dn) = —ds{6n)ds. 
The rescaling of the line element gives 6ds = K6nds, where k{z) is the local curvature of 
the boundary curve. At last, we have to vary the function / if it explicitly depends on the 
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k{z) = dn log 



contour. In particular, if this function is the harmonic extension of a contour-independent 
function on the plane, its variation on the boundary is given by 

= dr,{f{z) - f"{z))5n{z) , zedD 

This is an equivalent form of the Hadamard variational formula. (For more details see 
e.g. |2Sin2], where the Hadamard formula is extensively used.) 

The curvature. The local curvature of the boundary curve is defined as k = dd/ds, 
where 6 is the angle between the outward pointing normal vector to the curve and the 
X-axis. The formula 

w{z) 
w'{z) 

is an immediate consequence of the definition. For practical calculations, we also need 
the Laplace operator at boundary points in terms of normal (9„) and tangential (dg) 
derivatives, 

A = dl + d^,+ Kdn 

and the formula for the tangential derivative of the curvature, k' = dgK, through the 
boundary value of the Schwarzian derivative {w; z} = ^ — ( wT') • 

k'{z) = Im [u\z){w; z}) (8.2) 

{h> is the complex unit normal vector ()5.24|1 ). The variation of the curvature under small 
deformations of the contour is 

5k{z) = -{dl + K\z))5n{z) (8.3) 

Indeed, K + 5hi = where 59 = ds{5n) is the change of the slope of the normal 

vector. The line element is rescaled as ds —>■ ds + nSnds, so 6ds = n5nds. In the first 
order in 5n we then have k + 5k = n — d1{5n) — H?5n that is (j8.3|) . 



B. Variations of the trial functionals 

Let us consider in detail the most complicated case of the functional 

^^'^ = -^/ {(t>dn(t> + H)\dw\ 
27r J\w\=i 

To vary this functional, it is convenient to pass to the 2;-plane and introduce the function 

w\z) 



q{z) = log 



w{z) 



•4) 



It is harmonic in C \ D with the logarithmic singularity at infinity: q{z) = — log|2;| + 
0(l/|z|) as 2; — > oo. Obviously, q^{z) = log |w'(;z)| = q{z) + log 1^(^)1. In terms of this 
function, the curvature is k{z) = —dnq{z), and 

/(2) = — I (-qdnq + qe'')ds (8.5) 
2tc JdD 
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We would like to find the linear response of this quantity to small deformations of the 
contour. 

For clarity, let us deal with the two terms in (j8.5|l separately. We apply the rules 
given above and get, after some cancellations: 

S <j) qdnqds = <j> 5n \S/q\'^ds + <j>{5qdnq + qdnSq) ds 

where | Vgp = Using the Green theorem and the behaviour of the function 

q at infinity, one can see that the contributions of the two terms in the second integral 
are the same. The variation 6q at the boundary can be found by means of the Hadamard 
formula. The result is: 

6q = -e-'dnieUnf (8.6) 



Therefore, 



6 f qdnqds = ^ hvg|' - 2e''dn (e"'^a„g)'' 



6n ds 



where we used the Green formula again. The harmonic continuation is achieved by means 
of the identity e~'^dnq = 2TZe {^dzq) whose r.h.s. is explicitly harmonic and thus provides 
the desired harmonic continuation. Next, a straightforward calculation shows that 

eWn {e-^d^qf = -e'^d, (e-^sq) = {d^qf - dlq 

and we get 

5 j qdnqds = j (2dlq-\Wq\^ + 2k^) Snds 
The combination 2d^q — |Vgp can be transformed as follows: 

2d',q-\Vq\' = -2[,,\z)id',q-id,qf) + c. c 

where 2 {d1q—{dzqY) = {w; 2;} + | ^^'^ {"^S ^} is the Schwarzian derivative. Com- 

bining these formulas, we get: 

2dlq - I Vg|2 = - (z/2{w; z} + ly^{w]z} + ) (8.7) 
Now to the second term in ()8.5j) . We have: 

5 j> qe'^ ds = <j> Sndnqe'^ + <j> (e'^Sq + ge^ 5q) ds = 

= <j> 5ne^dn{q — q^) ds = — <j> 5n e^^ ds 

where we use the identity (9„log |w(-2)| = |w'(-2)| and take into account that \w'\ = e'^ at 
the boundary. Combining the results, we obtain: 

5/(2) = ^ ^ {i^\z){w; z] + z/2(z){w;z} - 2k^{z)) bn{z) ds (8.8) 

Variations of other trial functionals go in a similar way. Here we just list the results: 
b f \Vx\'d^z= i \Vx?5nds + 2 f VxV(5x)rf' 

Jd JdD Jd 
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^'z 



5 f Kxds= f {—dlx + ^9nX)^^ds+ f K6xds 

JdD JdD JdD 

S ^ xdnX"" ds = £^ (l VxP - {dnix - X""))') 6nds + 2 6x dnx"" ds 
5 / Axd'^z= f Ax6nds+ f dn{Sx)ds 

Jd JdD JdD 

To express the right hand sides through the variation of the potential W, one should 
plug the formulas ^J^, ^J^. 

C. The centrosymmetric potential at /3 = 1 

In the centrosymmetric case the potential W does not depend on the angular coordinate 
in the plane. We set W{z) = Wra.d{\z\'^) ■ At /3 = 1, the orthogonal polynomials technique 
is applicable. The symmetry of the potential implies that the orthogonal polynomials 
are simply monomials z"'. Therefore, the expression for the partition function simplifies 
considerably: 

N-l 

Zn = n 

n=0 

where 



h^= [ \z\^''e^^dh = 7c r x^'e^^^^^'^'^dx 
Jc Jo 



Introduce the function _ 

roc . 

h{t) = / e«(^-'*(")+"°§")cia; (8.9) 





then 

N-l 

logZN = N\og7i+Y, ^ogh{nh) (8.10) 

n=0 

We are going to find the ^-expansion of log Zj^ using the ^-expansion of h(t) obtained 
by the saddle point method and the Euler-MacLaurin formula. The necessary formulas 
are collected below. 



/i-Expansion of the free energy: a direct calculation. First of all we note that 
the density support D is an axially symmetric domain. We assume that it is a disk (not 
a ring) of radius R. The density of eigenvalues is a{z) = <7rad{\z\'^) ■ Clearly, 7rcrrad(a^) = 
—dx{xdxWrad{x)) . The radius i? is a function of to defined by the relation 

to = -R'WUR') (8.11) 

(here and below prime means f'{x) = df/dx). This relation follows directly from the 
definition 

to = / AWdh = -— (f dnWds 

4:71 Jd 47r JdD 

taking into account that dnW = 2RW^^^{R'^). 

The saddle point Xc in (j8.9p is determined by the equation 

t = -XcWU^,) (8.12) 
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Note the similarity with (jS.llI) : Xc is the squared radius of the disk for the fiUing n = t/h. 
We assume that there is only one saddle point. The standard technique (see (j8.16|) ) yields 
the following asymptotic expansion a.s h 0: 



hit) = ei(^^ad(x.)-H/;,,(x.)x.iogx.) / 2^xc ^ ^ ^^^2)^ ^g^^3) 
where pi is given by (j8.17|) with S{x) = PVrad(a;) + tlogx. It is easy to calculate: 

TTCTradfa^c) 



S"{Xc) 



Xr 



IV, ^^ -^(^'LA^c) , o^o^md(a;c) „vra,ad(a;c) 



Plugging this stuff into ()8.17|) . we get: 

In terms of the function x{^) = Xrad(|-2p) = | log(7r(Ji.ad) this expression acquires the 
form 



7rcrrad(2;c) V3 4 



_ 2 J_ 

q^radl^cj + 



The Euler-MacLaurin formula applied to (j8.1(jp yields 

logZ^v = ATlogTT + i r log/i(t)rft - ^log/i(to) + ^dt\ogh{t) 
n Jo 2 12 



+ C 

t=to 



where C is a constant which does not depend on Iq. (It can be found by a more detailed 
analysis around the point t = 0.) From ()8.13|) we have: 

log h{t) = ^{Wre.d{Xc) - W^/ad(a;c) Xc logxj + log ^271^ + 

+ ^logXc -XTUxc) + hpi + 0{h^) (8.14) 
The integral over t is transformed into the integral over Xc using dt/dxc = vrcrrad(2^c): 



/■to rJr! 

I logh{t)dt = 7r / aj-s_d{x) log h(t(x))dx 
Jo 



Finally, we obtain: 



f 

Fq = tx I (iyrad(a;) - W^/ad(3;) X \ogx) arsd{x)dx 
Jo 
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Fi/2 = - 77 / crrad(x) log(7rcrrad(a;))'ia; 
2 JO 

i^l = logi?^ - ixrad(i?^) - \r'}^...AR') + ^ ^(X:ad(^))'^^ (8-15) 

Reducing the 2D integrals in (j5.27p to ID integrals or values of Xrad, Xmd the 
boundary, 

Jd Jo 

f K-Xds = 27rXrad(-R^) 
JdD 

f Axdh = 4nR'xUR') 
Jd 

and taking into account that dnX^ = 0, we have: 

Fi = -1 logi?2 _ KUR') + xixUx)ydx + 2^iR\[^^{R') 

12 3 JO 

Comparing with (j8.15j) . we conclude that fi = —1/8. 



Asymptotic formulas. Evaluating the integral / en.^^^^dx by the saddle point method 
around the critical point Xc, 5"(xc) = 0, we get the asymptotic expansion as — 0: 



/ e^'^^Ux = e^^(-\/^^(l + hp, + h'p, + ...) 

J V \S"(Xr)\ 



S"ix,)\ 



where 



In particular, 



Pv 



2n 



27r(2n)! 



Pi 



S{x) — S{xc 



(X — Xr 



5iS"'{x,)f + 3\S"ix,)\S^''ix,)) 



2A\S"ix,)\^ 

In the text we also need the Euler-MacLaurin formula: 



Y: fin) = / fix)dx - -UXN) - /(O)) + -(/'(iV) - f (0)) + . . . 
n=o 2 12 



(8.16) 



(8.17) 



D. Derivation of the Polyakov- Alvarez formula 



In this appendix we outline the derivation of the formula ()4.11|) for the spectral determi- 
nant of the Laplace operator. We use the notation introduced in the main text. 

Consider a free Bose field X defined in a planar (compact) domain B and vanishing 
on its boundary. For simplicity we assume that B has topology of a disk. Let e^^dzdz be 
a background conformal metric on B. The classical action 



S = — I Xi-e-^'^A^Xe^^d^ 
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does not depend on the metric but the quantum theory does. The log of partition 
function of this field T = — |logdet (— c^^^Ab) represents the spectral determinant of 
the Laplace-Beltrami operator e'^-^As = 'ie~'^^^^^dzd2 acting on functions vanishing on 
the boundary (the Dirichlet b.c). Let z{w) be the univalent conformal map from the unit 
disk U onto B and = log \ dz/dw\. This map induces the conformal metric e^'^dwdw on 
U in the holomorphic coordinate w, with the conformal factor 

e2* = e2x|/(«;)|2 = e2(^+^) = V^ 

The Laplace-Beltrami operator in the coordinate w is 

4g-24.{^)^^^^_ The spectral deter- 
minant of this operator can be expressed in terms of $ and the curvature of the boundary. 
The result is referred to as the Polyakov- Alvarez formula [211 EO] • Below we give a short 
derivation of this formula. 

Variation of the partition function over the metric introduces the trace Tzz = {\dX\^) 
and the holomorphic component T = (((9X)^) of the stress energy tensor: 

The trace of the s.e.t. has bulk and boundary parts, and so has JF: 

Tz-zdzdz = T'^f'dzdz + T^^J^^^'^^'^yrfs, ^ = ^^"l"^ + ^boundary ^g^Q) 

The holomorphic component is continuous across the boundary in a regular fashion. The 
components T^z and T are related by the conservation law 

d{y/gTz-z) + y/gdT = (in the bulk) (8.21) 

2dsT^t'''^'''y + Im{u^T) = (on the boundary) (8.22) 

where u is the unit normal vector to the boundary. 

The strategy to compute JF is as follows. First we compute the holomorphic compo- 
nent T, then find Tzz through the conservation laws (|8.21|) . (|8.22|) . Finally we integrate 
eq. (IHini). 

The simplest way to proceed is as follows. Since the holomorphic component is 
continuous across the boundary, it depends only on the overall metric y/g = e^-^"*"^*^. 
Therefore, it is sufficient to compute T at x = 0. In the coordinate w, the domain is 
the unit disk, so T{w){dwY vanishes. The conformal transformation w — > z generates 
the metric e~'^^ and transforms the holomorphic component of the stress energy tensor 
as T — > T + -^{w; z}. Therefore'', 



1 1 

T(z) = — \w; z} = 

at X = 0, or, more generally, 

1 



(d\og\dz/dwf^ + 2d'^ log \dz/dw\ 



T = ((Slog + 2d' log ^) (8.23) 



^One can obtain this textbook result by the direct computation of {dX{z)dX{z')) cx dzdz'G{z, z') 
and extracting the finite part of the resuh as z' — > z with the help of the explicit formula for the Green 
function KT^ . 
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at X 7^ 0. The next step is to find the trace of the stress energy tensor. The bulk part 

7^'"^' = ^Alogv/^=^Ax (8.24) 

follows from 1)8.211) and represents the gravitational anomaly. It is proportional to the 
scalar curvature of the metric in the bulk and does not depend on the shape of the 
domain. This gives (after integrating eq. ()8.19)l ) the textbook result for the bulk part of 
the spectral determinant: 

= ~^/u^°Sv^Alogy^d2«; (8.25) 

We will use the same trick in order to compute the boundary part. Let us first find it 
at X = 0. Comparing ()8.22|) and ()8.2|) we conclude that it is ^dgH, where k is the local 
curvature of the boundary. This gives the boundary part of the gravitational anomaly: 
= —^i^- At X 7^ one should simply add the boundary curvature generated 
by the metric e^^ to the local curvature k: 

T^r'^^'^yds = -^i^ + dnX)ds . (8.26) 

Using the formula hzds = {dn4>+e~'^)ds = dn(pds + \dw\ (equivalent to (HHl)), the boundary 
contribution to the gravitational anomaly can be written as 

^boundary = - 1 {2\dw\ + \og ^ ds) 

This form is ready for integration with the result 

^boundary^ 1 / 1 J | | + ^ i ^ \dw\ + i d^log ^ \dw\ 

967r Jdu 24:71 Jdu Wtt Jd\j 

(8.27) 

where we take into account that dn on dD is e~'^dn on the unit circle. Combining with 
fl8.25|) . we obtain the finite part of ()4.11|) (recall that log^^ = 2$). 
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